Multispin correlations and pseudo-thermalization of the transient density matrix in 
solid-state NMR: free induction decay and magic echo 
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Quantum unitary evolution typically leads to thermalization of generic interacting many-body 
systems. There are very few known general methods for reversing this process, and we focus on the 
magic echo, a radio-frequency pulse sequence known to approximately "rewind" the time evolution 
of dipolar coupled homonuclear spin systems in a large magnetic field. By combining analytic, 
numerical, and experimental results we systematically investigate factors leading to the degradation 
of magic echoes, as observed in reduced revival of mean transverse magnetization. Going beyond the 
conventional analysis based on mean magnetization we use a phase encoding technique to measure 
the growth of spin correlations in the density matrix at different points in time following magic 
echoes of varied durations and compare the results to those obtained during a free induction decay 
(FID). While considerable differences are documented at short times, the long-time behavior of 
the density matrix appears to be remarkably universal among the types of initial states considered 
- simple low order multispin correlations are observed to decay exponentially at the same rate, 
seeding the onset of increasingly complex high order correlations. This manifestly athermal process 
is constrained by conservation of the second moment of the spectrum of the density matrix and 
proceeds indefinitely, assuming unitary dynamics. 
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I. INTRODUCTION 

Nuclei in solids often enjoy relatively weak coupling 
to the environmental, non-magnetic degrees of freedom, 
marked by exceedingly long thermalization times, e.g. 
T\ that can reach several weeks^. Under these condi- 
tions coherent quantum many-body dynamics whose in- 
trinsic time scales are typically some microseconds can 
be studied and manipulated with relative ease. Much 
of the progress in elucidating local physics and chem- 
istry of diverse substances over the past six decades has 
been achieved by matching detailed solutions of few-body 
dynamics with controlled experimental studies whereby 
through a combination of DC and RF fields desired 
many-body evolution can be realized to a high degree 
of precision. The generation of so called "echoes' — is a 
particularly central tool in the NMR arsenal. All echo 
schemes strive to pattern the finite time quantum evolu- 
tion operator in such a way as to render it equal to the 
identity operator, thereby returning the many-body sys- 
tem to its initial state. Most famously, the Hahn echo^ 
recovers the transverse magnetization dephased by local 
field inhomogeneities (and chemical shifts) with a sim- 
ple reflection about an axis in the plane of spin pre- 
cession. Remarkably, such single-spin corrective time- 
reversal action can be generalized with considerable suc- 
cess even to cases where the effective chemical shifts are 
time dependent (but mutually uncorrelated)^. Not sur- 
prisingly, time-reversing correlated multispin dynamics 
is more subtle as the effective Hilbert space is no longer 
bounded. To our knowledge there is no general system- 
atic protocol for achieving this - along with Maxwell's 
demon such general Loschmidt echoes^ were believed 



to be purely of conceptual significance. Yet, for dipolar 
spins subjected to strong DC fields an explicit solution 
to this problem, the so called "magic echo," was found^, 
and has since served as the foundation for NMR stud- 
ies of correlated spin motion, including measurements of 
imaging^, spectroscopy^, spin-diffusio n 10 ! 11 , to investi- 
gate nonlinear dynamics in highly magnetized liquid*^, 
and, as in this work, for phase encoding and observing 
multispin correlations in non-commuting basest 3 -. 

Our work pursues two different lines of inquiry which 
lead to a complementary characterization of the magic 
echo protocol. First, the experimentally observed efficacy 
of the magic echo, i.e. the relative amplitudes of the net 
refocused magnetization acquired using different magic 
echo conditions, are compared with average Hamiltonian 
theory extended to second order, and also to results of 
numerical simulations. Here, we specifically focus on the 
relative importance of the finite width of the RF pulses 
vs. the total RF power applied". Second, the time 
dependence of the full density matrix is characterized 
by measuring the onset and spreading of multispin cor- 
relations (also known as "spin-counting" spectroscopy) 
both during the usual free induction decay (FID, i.e. af- 
ter a 7r/2 pulse from equilibrium) and following magic 
echoes with a few different, relatively long refocusing 
times. Interestingly, while these initial states are dis- 
tinct and therefore evolve differently at short times their 
late time evolutions are apparently very similar, in the 
sense that certain rates of growth and decay are quantita- 
tively indistinguishable. In the balance of the Introduc- 
tion we summarize the basic physics of the magic echo 
and spin counting technique to help an impatient reader 
skip the following two Sections where a detailed discus- 
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sion of notation and the theoretical framework (Section 
II), as well as experimental and numerical methodology 
(Section III), are given. Results are presented and com- 
pared with theoretical expectations in Section IV with 
conclusions and discussion of future directions of research 
relegated to Section V. 

The magic echo involves two key ingredients. First, the 
application of a strong DC field, B , along the z direction 
simplifies the interaction between two nuclei to 

Hd = £jr (1 " 3cos 2 9) (2I°I° - \{ip k + Ijl£)) 

= D jh [2l9%- (1) 

where 9 is the angle between the internuclear vector rjk 
and the applied magnetic field which is oriented along 
the z-axis 1 ^, 7 is the gyromagnetic ratio of a nucleus, 
and J ,/* are canonical spin operators. To set the no- 
tation for the remainder of the discussion the single spin 
(Larmor) precession frequency corresponding to the DC 
field is cjl = 7 -Bo, while the internal dipolar fields can 
be associated with "dipolar frequency" wd = "f 2 /r 3 (r 
refers to the smallest internuclear distance). Typically 
ojd ~ 10 4 Hz, while wl ~ 10 s Hz so the rapidly oscillat- 
ing terms omitted from Eq. 1 arc effectively suppressed 
by a small factor ojd/wl- Also, in what follows we will 
always work in the frame rotating at wl (with respect 
to the lab frame). Second, continuous application of a 
strong resonant field (i.e. at Hi, in the x — y plane 
"locks" the magnetization as the spins execute approxi- 
mate Rabi oscillations at a frequency u>i = 7 Hi. In this 
Rabi frame of reference, co-rotating with the RF field, 
the dipolar interaction can be further reduced, approxi- 
mately (i.e. assuming ujd/uji -» 0), to 

-^[27^-l(J;4 + /7J-)], (2) 

where J ,/* refer to the quantization axis in the sec- 
ond rotating framed, i.e. if Hi is along x then 1° = I x . 
Most importantly, the relative minus sign of the (0,0) 
and (±, =f) terms between Eqs. 1 and 2 translates into an 
overall sign reversal of spin interactions in the presence of 
the resonant field. After a few relatively minor additional 
steps the entire propagator during the application of the 
RF field (often referred to as an "RF burst" ) becomes 

Ui = e+ m D t/(.2 h ) t (3) 

thus allowing the effective reversal (also referred to as 
"refocusing" or "time-suspension") of the entire many- 
body trajectory, which will typically consist of evolution 
under Uu = e~' inDt / h as well as U\. 

The idealized discussion above relies on a number of 
assumptions, of which the most tenuous is the availabil- 
ity of a very strong and precisely applied RF field. Quite 
generally, the efficacy of magic echoes is expected to de- 
grade when longer trajectories are reversed either due 



to corrections beyond the zeroth order approximation in 
wd/wi (in the so-called average Hamiltonian theory )^2r— 
or pulse artifacts (e.g. ring-down effects). One gener- 
ally expects the former (latter) limitation to dominate 
at small (large) RF field strength, albeit for very dif- 
ferent reasons (e.g. pulse artifacts are of purely experi- 
mental origin). Surprisingly, very little is known about 
the degradation of magic echoes in the limit of long time 
suspension, e.g. > 1 ms, and/or intermediate RF fields. 
Although it is known how to achieve even longer time 
suspensions using short multiple (hundreds or even thou- 
sands) magic echoes in sequence^^, some applications, 
e.g. the study of spin diffusio n 10 ' 11 , can benefit from the 
availability of a long duration single magic echo. Un- 
derstanding and eventually negating the factors leading 
to the degradation of long magic echoes, especially under 
realistic conditions of fixed total RF power applied to the 
sample, is a long term goal. 

While revival of the readily observable (uniform) mag- 
netization is customarily taken as a proxy for the quality 
of the echo, we expect, at least in principle, the entire 
density matrix to be returned to near its initial (a 7r/2- 
rotated thermal) value p(t = 0) ~ 1 - Y,j If ■ I n the 
remainder of this work, we omit the identity term from 
p(t = 0), as is customary in NMR at infinite tempera- 
ture. What can be said about the structure of the density 
matrix? The density matrix of nuclear spins in a rigid 
lattice during a conventional FID is described by the evo- 
lution of single-spin, single-quantum coherence (i.e. £ If 
term) into single-quantum, multiple spin correlations (i.e. 
terms involving IfI^...If) . As these multispin correla- 
tions are not readily observable, the density matrix may 
appear to have thermalizcd but this is forbidden, strictly 
speaking, by the exact conservation of its spectrum un- 
der unitary evolution. As we will explain below, the ex- 
periment used to measure multispin correlations in fact 
accesses the second moment of the spectrum of the den- 
sity matrix (the so called purity) and details the process 
of pseudo-thermalization, by which athermal correlations 
propagate (or spread) away from simple and readily ob- 
servable parts of the density matrix to complex higher 
order correlations. 

The early onset of these multispin correlations in a 
crystal has been measured in recent experiments^ using 
a previously developed phase encoding technique. The 
basic trick essentially amounts to introducing steps to 
rotate the spin axis by a controlled amount thereby con- 
verting single-quantum operators into a linear superposi- 
tion of operators of different coherence order, which can 
be distinguished. We extend this sort of study of the 
"anatomy" of the density matrix both to the late time 
regime for the usual FID but also, more interestingly, to 
post magic echo evolution. Direct comparison of short 
time dynamics reveals quantitative differences between 
the two types of the density matrices already at short 
times. Interestingly, and perhaps surprisingly, long time 
dynamics appears to be more universal. In addition to 
the "sum rule" mentioned above that guarantees that 
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FIG. 1. The magic echo sequence used in this manuscript. 
After an initial 7r/2a, pulse, the spins evolve for a time r, fol- 
lowed by a sandwich consisting of two RF bursts along x and 
x, each of duration 2t between two 7r/2 pulses along y and 
y. Because of the timing of this sequence, the echo appears 
at a time r after the final 7r/2 pulse. For the average Hamil- 
tonian calculations the initial 7r/2 x pulse is not considered. 
The toggling frame Hamiltonians during time periods A-F 
are provided in Tabled] 



the growth of higher order correlations must come at the 
expense of lower order ones, we observe a high degree 
of similarity in the dynamics. For a sufficiently short- 
ranged Hamiltonian (which appears to include dipolar 
interations^) both of these processes (decay and spread) 
appear to be simple exponential in time, and we docu- 
ment the existence and correspondence of two multiple 
spin correlation decay rates in the same experiment and 
across different initial states^. 



II. THEORY OF THE MAGIC ECHO AND 
SPIN-COUNTING SPECTROSCOPY 

In this Section we review and extend the average 
Hamiltonian theory of the magic echo, isolating finite 
pulse width effects from other ( "burst" ) corrections con- 
tributing to an echo's demise. We also derive explicit 
expressions connecting the phase encoding experiment 
with density matrix elements and its purity. All of this is 
available elsewhere, except, to the best of our knowledge, 

I 
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where r c = 6r + 2i 7I y 2 f° r the magic echo (see Fig. [TJ. 
Further analysis is facilitated by defining the so-called 
"toggling frame Hamiltonian" which in the absence of an 
RF field is rotating at the Larmor frequency ujl about z 
and is otherwise rotating about the RF field at u\, with 
r » t„i2 = ^j-, the duration of the ir/2 pulses. We list 



the first and second order average Hamiltonian calcula- 
tions for the magic echo and the extension of the spin 
counting experiment to probe multiple spin correlations 
following a magic echo. 



A. The Magic Echo 

In a large, static magnetic field the spin dynamics in 
solid state NMR of spin- 1/2 nuclei are dominated by the 
secular dipolar Hamiltonian T-Ld ■ To describe the effect of 
time dependent fields, average Hamiltonian theory may 
be usecU^. In this formalism, the evolution of the density 
matrix is 

e (r c ) = UdURFUdQiOKUdURpUd)- 1 

^u{ Tc ) e {0)u-\T c ). (4) 

The propagators Urf and Ud refer to evolution with and 
without the RF field, respectively, and r c is the time for 
a complete magic echo cycle. The propagator U rf may 
be written by the Dyson series 

U RF (t)=Texp(-i f U RF {ti)dt{). (5) 
Jo 

In the above expression T is the time ordering operator 
and % rf is the RF Hamiltonian in the rotating reference 
frame (which is explicitly time-dependent). The object 
of average Hamiltonian theory is to impart a well-defined 
time dependence on the internal interaction, which is oth- 
erwise time independent, via a sequence of experimen- 
tally controlled RF pulses. Over a cycle of RF pulses of 
duration r c , the propagator U may be written in terms 
of the Magnus expansion^ 

C/(r c ) = e-^° + « I+ -K (6) 

The zeroth, first, and second order terms in the Magnus 
expansion are 

(7) 
(8) 
(9) 

I 

the toggling frame Hamiltonians % for the magic echo 
sequence shown in Fig. [TJ in Table |TJ 

In our analysis of the average Hamiltonian we only 
show interactions between pairs of spins. Because of the 
commutators in Eqs. [S] and El operators involving three 
or more spins generally appear in the higher-order tcrms- 
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= ~ f TC dt 3 [ t3 dt 2 r t2 ([H(i 3 ),[H(i2),H(ii)]] + [H(ii),[H(t 2 ),H(i3)]])^i, 
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TABLE I. The toggling frame Hamiltonians for the magic echo pulse sequence depicted in Fig. [T] Here 9\ = 2twi is the flip 
angle of the first x RF burst as seen in Fig. [T](time period C). 

Time period Duration Toggling frame Hamiltonian 



A 
B 

c 

D 

E 
F 



£>i2(3/r/!-/i-i 2 ) 

ZD 12 [im cos 2 (wit) + (7f7£ + If 1$) cob (wit) sin (wit)+ 7f 7f sin 2 (wit)] - D 12 Ji • / 2 
3£>i 2 [7f 7f cos 2 (wit) - (7f 7| + 7j7f ) cos (wit) sin (wit)+ 7»7f sin 2 (wit)] - TWi ■ 7 2 
373i 2 [7f7f cos 2 (-wit + 6»i) + (7f7| + 7~7f ) sin (-wit + 9i) cos (-wii + 0i)+ 7f7| sin 2 (-wit + 6>i)] 
-Di 2 7i -7 2 

tt/(2wi) 3Di2[7f7| sin 2 (wit) + (7?7f + If 7|) sin (wit) cos (wii)+ 7f7f cos 2 (wit)] - Di 2 7i ■ 7 2 
Di 2 (3/ 1 z 7|-7 1 -7 2 ) 



T 

K/2 = 1"/( 2w l) 

2r 
2r 



lTr/2 = 



we neglect these terms for simplicity. The first three 
terms in the Magnus expansion for the secular dipolar 
Hamiltonian in the magic sandwich shown in Fig. [T] are 

I 



given by (note that Ref. [l8| does this zeroth order average 
Hamiltonian calculation for a slightly different variation 
of the magic echo sequence) : 



1/0 

H D 



-Dl2*7r/2 



I 1 -I 2 -3^7| + -(Jf7f +7f/|) 



3D 



12 



sin 2 (2wir) 



(10) 



9L> 2 2 sin 2 (2rwi) sin (4rwi)(/ 1 2 + 7|) 
16w 2 r c 



(11) 



- 2 9£>i 2 t W2 £> 2 sin(4Ta;i) 
«rj = L^-(drwi + 

l07TT c 2 



+ ^|(__L +T 2.= 
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1 

16 

Tg 



- 3toji 



cos(4twi) 
8 

7rcos(4rwi) 



)(47f7f-27f7f-27?7j) 
cos(8toji) 



16 

7rcos(8ra;i) 
16 



+ twi sin(4rcji))(47j c 7| 
rwi sin(4rcLii) 



)(27 r 7*-27f7f)]. 



(12) 



r 



In the limit of infinite wi the above expressions vanish, as 
expected. With 71^=0 (n = 0, 1, 2 . . .), one refocuses the 
many-body spin dynamics resulting in the observation of 
a magic echo. The zeroth order term term contains im- 
perfections from both pulse width effects and finite RF 
during the burst. The first order term, Hp, results from 
the commutator between the toggling frame Hamiltoni- 
ans during times C and D in Fig. [1] and is, therefore, 
free of pulse width effects. This term vanishes if the du- 
rations of the RF bursts are such that 2_wit is a multiple 
of 7r and/or if D12/UJ1 -> 0. Even if + 0, it cannot 
degrade the magnitude of transverse magnetization - it 
involves only 7^ z operators and only produces a rotation 
about the z-axis. The second order term in the Magnus 
expansion, "H^, arises only from commutators including 
time periods B and/or E (the 7r/2 pulses of the magic 
sandwich) - this term vanishes for infinitely narrow (5- 
function) 7r/2 pulses in the magic sandwich. Importantly, 
for finite t„/2 it includes terms which increase with echo 
time t. Thus, to second order finite RF pulses will in- 
duce imperfections in the refocusing offered by the magic 
sandwich with increasing r which may be evident in both 
the multiple spin correlations and in the single-quantum 



signal from a magic echo. 



B. Multiple Spin Correlations 

When NMR signals are detected by inductive coupling 
to a coil the measured signal at time t is given by 



S(t)~Tr[p(t)7 + ]. 



(13) 



Only single spin, single quantum terms survive the trace. 
The growth of multiple spin single quantum coherence 
can be observed in the short time expansion of the density 
matrix 

g{t) = g{Q) + U [g(Q),B]~ [[g(0),H] , 77] + . . . , (14) 



where 



e(o)~-£/* 



(15) 



after a 7r/2 y pulse (i.e. <j) = in Fig. [2]). Using Eq. [I] for 
the Hamiltonian the evolution leads to 
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with i within the time interval labelled T in Fig. [2j The 
insight into how to begin resolving the different terms 

I 



It is clear that upon basis rotation single quantum and 
zero quantum operators are rotated into each other, pro- 
ducing operators changing polarization along x by ar- 
bitrary amounts. The growth of such observable x- 
multiquantum operators can be taken as a proxy for 
the onset of multispin correlations (in fact, the two are 
uniquely related). All that is left is to "count" coherence 
orders in the non-commuting basis, which is done by in- 
troducing phase-modulated rotations, as in Fig. 2 and 
Eq. 19. In the experiment shown in Fig. [2] the transfor- 
mation of basis is performed using the first, fourth, and 
fifth 7r/2 pulses. To encode odd orders, the first pulse is 
phase shifted by 7r/2 from what is shown in the figure. In 
the present work, with the magnetization in the x basis, 
the phase cf> is varied from to Air with 120 total values 
of <j> to encode up to the ±30th order of spin correlations. 
The zeroth order spin correlations have coherence order 

I 



where t" is the time after the last pulse and H is the 
dipolar Hamiltonian. This propagator explicitly assumes 
that the second magic sandwich perfectly refocuses the 
dynamics (this assumption is valid as this magic echo is 
relatively short). The first magic echo will be absent if 
one is studying multispin correlations in the FID. With 
t' = T' - 1" - S the peak of the echo is located at t' = T 
while the propagator simplifies to: 



[U+]- 1 = < " :J < ""' (19) 



(16) 

I 

in the density matrix comes from simply rewriting it in 
another, non-commuting basis, e.g. polarized along x 



(17) 

I 

in the a;-basis, corresponding to the first term plus four 
of the six terms in the quadratic-time-dependent portion 
of Eq. \T7\ as well as those of coherence order zero from 
higher orders in time. The higher-order terms are con- 
verted back into single-spin terms by a magic echo, i.e. 
the final magic echo of Fig. [2] For signal processing, 
one implements a Fourier transform of the resulting sig- 
nal (with respect to <f) to generate a spectrogram of the 
coherence orders. Because of the mixing of coherence 
order N among spin correlations with N or more spins, 
the signals do not correspond directly with iV-spin single- 
quantum coherence terms from Eq. [16] but instead rep- 
resent the correlations among at least N spins. Defining 
a shorthand for the total moment I = Y,j Ij we can write 
the (inverse) propagator of the pulse sequence in Fig. [5] 
explicitly 



(18) 

I 

and the observed signal is 

s* = (u; 1 i'U4i x ) 

I iHt' -i<pr -iHT jx iHT id>r -iHt' T x\ 

= (e e v e I e e v e 1 ) 

= {m\e- i + I '' 2 I x (T)e i + I '/*\ n ) (n\e^ r ' 2 I x (t')e-^ r ' 2 \m) 

= E e i{n - m) *[I X (T)] mn [I x {t')] nm E^„e'"* (20) 

n,m n 

Note that in addition to assuming a perfect magic echo 
this protocol also assumes perfect conservation of I z , 
which is valid in the high field limit. Clearly, £ n S n = 
Ti[p 2 } for t' = T. 
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FIG. 2. Pulse sequence used to observe multiple spin correlations following a magic echo. For the observation of multiple 
spin correlations following the FID, the first magic sandwich is removed as are the delays r before and after it. The time 
suspension sequence is used to suppress artifacts and suspend evolution under the dipolar interaction^. The delays T' and 5 
are adjusted so that the refocusing of the dipolar interaction peaks after instrumental ringdown following the end of the final 
magic sandwich. The flip angle 9y = oj\T' . 



III. METHODOLOGY 

A. Experimental procedures and sample 
characterization 

Experiments were performed in adamantanc and cal- 
cium fluoride at room temperature. Adamantane is a 
plastic crystal in which molecular tumbling averages out 
intramolecular interactions 2 ^, leaving every proton to in- 
teract with each of the 16 protons in all of the 12 nearest 
neighbor molecules in addition to those further away. Be- 
cause of tumbling these molecules act as if all the protons 
in a molecule are at the center of that molecule. The ef- 
fective dipolar coupling D c h in adamantane is « 4.7 x 10 
rads/s 2 ^, and is determined using a measurement of the 
linewidth. The calcium fluoride sample used in this study 
was a single crystal oriented such that the duration be- 
tween the first two zero crossings in the FID (note that in 
Figs. [3] and |4] only the magnitude of the signal is plotted) 
is approximately 40 /is, which lies between the values for 
the [110] and [111] axes along the Zceman field. 1 H and 
19 F NMR experiments were performed at 179.445 and 
168.824 MHz using a Tecmag Apollo spectrometer. 

In Figs. [3] and [4] we show a series of magic echoes and 
FIDs in adamantane and calcium fluoride respectively. 
The data are normalized and time-shifted to emphasize 
the overlap in the long-time portion of the decays 2 ^ In 
the case of the magic echoes, the initial portions of the 
signals vary with echo time but eventually decay and 
oscillate at the same rate. We show the results g and 
fi of a fit of the long-time portion of the decays to 

F(t) = F exp (-gt) sin (fit + $) (21) 

in Table |TT] for adamantane and calcium fluoride. This 
form for the long-time portion of the FID has a long 
history^ and has recently been predicted using the no- 
tion of microscopic chao s 21 ' 27 . The sinusoidally mod- 
ulated exponential characteristic has also been ob- 
served for the FID when measured out to six orders of 
magnitude 2 ^. We note that our measurements for cal- 
cium fluoride do not correspond to those in Ref. H^, be- 
cause our crystal is not exactly in the [110] orientation. 
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FIG. 3. (Color online) Experimental FID and magic echoes 
measured in adamantane. To overlay the data, we have time- 
shifted and normalized each echo decay. Even though the 
initial decays vary significantly, they all approach the same 
long-time decays. In this figure the points represented are: 
the FID ("/black) and magic echoes at different t,t— 100 /is 
(♦/red), r=150 fis (»/blue), and t= 200 /is (y/green). 



Figure [5] highlights a representative experimental spec- 
trogram for spin correlations measured in adamantanc 
with a magic echo time r = 100 fia and an evolution time 
T = 240 fis. In all experiments performed, we used 4 
signal averages with a recycle delay of 5 (30) seconds 
in adamantane (CaF2). In our experiments, S was set 
to 10 /us and T = T + 90(70) n& in adamantane (CaF 2 ) 
to avoid probe ring down artifacts. Our experimental 
pulse sequence also implemented a simple phase cycle to 
remove baseline artifacts and imbalance in the receiver 
channels 2 ^. The time of the echo peak for the = 
phase encoding step was used when computing the spec- 
trogram. 
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TABLE II. Values of the parameters f2 and g determined for the fit of the FID and magic echo signals to the equation 
Fq exp (~gt) sin (Sit + <3>) in adamantane and calcium fluoride. The parameters shown in the tables (and associated error bars) 
were determined from six independent fits to the data. 



adamantane 


SI (rads/ms) 


g (ms *) 


CaF 2 


SI (rads/ms) 


g (ms 1 ) 


FID 


30 ±1 


32 ±4 


FID 


71 ± 1 


34 ±3 


100 fis magic echo 


31 ± 1 


31 ±3 


50 /is magic echo 


73 ±2 


37 ±2 


150 [is magic echo 


36 ± 1 


33 ±2 


100 fjs magic echo 


74 ±2 


33 ±3 


200 [is magic echo 


33 ± 1 


34 ±3 


130 fis magic echo 


75 ±3 


33 ±3 



10 n 




Time (|is) Coherence Order 



FIG. 4. (Color online) Experimental FID and magic echoes 
measured in calcium fluoride. To overlay the data, we have 
time-shifted and normalized each echo decay. Even though 
the initial decays vary significantly, they all approach the 
same long-time decays. In this figure the points repre- 
sented are: the FID ("/black) and magic echoes at differ- 
ent r, r=50 /is (*/red), r= 100 /is (*/blue), and r=130 /is 
(y/green). 



B. Notes on simulations 



FIG. 5. Representative even order experimental multiple spin 
correlations, obtained in adamanatane for magic echo time 
r = 100 fj,s and evolution time T = 240 /xs (see Fig. [2] for 
definitions of r and T). We performed 4 signal averages and 
incremented <j> from to 4-zr by 7r/30 (120 phases <j> were used). 
The spectrogram shown in the figure is obtained by a Fourier 
transform of the echo peak amplitude variation with phase. 
The growth and subsequent decay of coherence peak ampli- 
tudes for adamantane and calcium fluoride are plotted in Figs. 
E]and H 



Numerical simulations were performed by integrating 
equations of motion for the expectation values of spins 
in the classical limit, with quantum spin operators re- 
placed by "classical spins," i.e. three component unit 
vectors obeying the appropriate Poisson bracket. These 
deterministic Bloch equations were integrated using a 4th 
order Runge-Kutta method with a fixed time step, At. 
Optimally, this time step is chosen to be short enough 
to sample single spin precession at least ~ 10 times per 
typical precession period but not too short to overbur- 
den the computation. Going over and performing the 
entire simulation in the rotating frame (which is justi- 
fied in the large DC field limit) the value of At is set by 
the strongest, nearest neighbor, dipolar coupling, unless 
the RF field is on, at which point At is reduced accord- 
ingly, as the precession about the RF field is faster than 
in dipolar fields. To facilitate simulations we have trun- 
cated the dipolar Hamiltonian to retain only the nearest 



neighbor couplings and use cubes with 16 3 spins. Of the 
various possible finite-size, -range, and -time effects we 
have found empirically that finite time effects controlled 
by At are by far the most serious and focused on minimiz- 
ing them. We spot-checked that quoted results are robust 
against increasing both the interaction range and lattice 
size. Lastly, to overcome the low signal-to-noise ratio in 
such small lattices we artificially reduce the temperature 
so that j3 = (ksT)' 1 = 0.2 (and we average over 100 ini- 
tial conditions) . Since this value of magnetization is five 
orders of magnitude higher than that of CaF2 at room 
temperature in high fields, we explicitly checked that the 
spin dynamics is still effectively governed by infinite tem- 
perature correlations by working at higher temperature 
(and thus at lower mean magnetization) and only observ- 
ing enhancement of statistical noise. 
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IV. RESULTS AND DISCUSSION 

A. Degradation of the magic echo 

In Figs. [6]and[7]we show intensities of the magic echoes 
in experiments using adamantane and in simulations, re- 
spectively, versus the RF field strength power applied 
during the 7r/2 pulses uii/D e s (D e g rs ID 12 in a cubic 
crystal). In the experiments, whose results are shown in 
Fig. [5J the amplitude of the initial 7r/2 and all RF bursts 
was fixed to 8.1D e g. The experiments clearly show a 
decrease in the echo amplitude with both echo time r 
and with decreasing RF amplitude u>i of the ir/2 magic 
sandwich pulses. In Fig. [7] magic echo amplitudes for 
varying RF amplitude and r arc shown for simulations 
of classical spins. In these simulations the amplitudes of 
the RF bursts (during periods C and D) are kept con- 
stant at 6.5-Deff rads/time. The data reveal a decrease in 
the echo amplitude with increasing r but the dependence 
on cjj is weaker than that observed in the experiments. 
These classical simulations became unstable for longer 
relative echo durations so the quantity tD c k was larger in 
the experiments than in simulations for the largest echo 
times we used. When TD c g = 5.7 the variation in the 
magic echo amplitude decreases by ~ 21% when u}i/D c g 
is reduced from 6.5 to 3.3. This compares to the exper- 
imental case where TD e g = 4.7 and the echo amplitude 
decreases by ~ 4% when Wi/-D e ff is reduced from 6.6 to 
2.7. Coupled with the weaker relative RF in the simu- 
lations' spin locking periods the classical simulations at 
least qualitatively capture the features observed in the 
experiments on adamantane emphasizing the contribu- 
tion of finite pulse widths in the degradation of the magic 
echo cycle. Most importantly, the relative influence of fi- 
nite width 7r/2 pulses becomes diminished at any r with 
only moderate decrease of t^n- In summary, we ex- 
perimentally observed a degradation in the magic echo 
efficacy with decreasing RF amplitude of the magic sand- 
wich 7r/2 pulses which becomes pronounced for large echo 
times. There is also a flattcning-out of the magic echo 
refocusing for large RF amplitude as increasing the RF 
amplitude does not provide better refocusing. While this 
observation may be due to instrumental artifacts such as 
phase transients (which are known to increase with in- 
creasing RF amplitude), the fact that it is also borne out 
in numerical simulations suggests that this flattening is 
an intrinsic effect, in principle captured by the average 
Hamiltonian treatment. 



B. Growth, decay, and spread of multispin 
correlations 

In Figs. H] and Owe show the various multiple spin cor- 
relations plotted as a function of the echo spacing r as 
well as the evolution time T observed in the experiments. 
Each echo time's multiple spin correlations in Figs. [5]and 
[5] have been renormalized separately (one rcnormaliza- 
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FIG. 6. (Color online) Intensity of the magic echo observed in 
adamantane for different values of r and oji corresponding to 
the RF field applied during the 7r/2 pulses of the magic sand- 
wich. As described in the text, the RF amplitude during the 
spin-locking pulses (periods C and D in Fig. [1} was kept con- 
stant for all experiments. The values of r are 100 /is ("/black), 
200 [is (*/red), 250 /is (*/blue), and 350 ^s (y/green) corre- 
sponding to tD c h = 4.7,9.4,12, and 17 respectively. As dis- 
cussed in the text, the variation in the echo amplitude with 
loi and r is due largely to the second and higher order terms 
of the Magnus expansion. 



tion per echo time r, for all orders of spin correlations). 
The adamantane FID data are in good agreement with 
a previously published study^. The agreement between 
the CaF2 data and previously published data^ is rela- 
tively close for short times, although should be noted that 
in the authors of Ref. [TH normalize the coherence orders 
so that at every point in time the sum of all orders is 
one (the correspondence is apparent if one plots the data 
this way). However, there are some differences, which be- 
come especially apparent in the higher order correlations 
at longer times. The orientation of our crystal appears to 
be different when the FIDs are comparcd-prior data was 
acquired with the crystal's [110] axis along the magnetic 
fields. The initial growth of the multiple spin correla- 
tions details a significant difference between the FIDs and 
the various echoes, revealing imperfections in the time re- 
versal of the magic echo. We now turn to a comparison 
of the FID and magic echo data. Quite generally, higher 
order correlations appear large in magic echo traces as 
compared to those of the FID, which is not surprising. It 
is also apparent that the higher order correlations peak 
earlier for the echoes than for the FID, an additional in- 
dication of imperfect refocusing and leftover multispin 
correlations at the peaks of the magic echoes which are 
not present in the FID. 

Looking beyond the initial growth stage, Figures[5]and 
IHl reveal that the multiple spin correlation amplitudes dc- 
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FIG. 7. (Color online) Magic echo amplitudes for lattices of 
classical spins as a function of wi during the n/2 pulses of 
the magic sandwich. As in Fig. [6] the RF amplitude during 
the spin-locking pulses (periods C and D in Fig. [1]) was kept 
constant. The values of r are 0.2 ("/black), 0.5 (*/red), 1 
(*/blue), 2 (v/green), 4 (< /orange), and 8 (►/purple) corre- 
sponding rD e s = 0.14, 0.36, 0.71, 1.4, 2.8, and 5.7, respectively. 
The effective dipolar coupling strength is D c ff = 0.71, as de- 
termined by the linewidth. Note that the values of rD e g are 
smaller here than nearly all the values in Fig. [5] and that, 
correspondingly, the variation in echo amplitude is smaller. 
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FIG. 8. (Color online) Measured even-order multiple spin 
correlations in adamantane as a function of time T in Fig. 
[2] The shapes represent the FID (■) and echoes at various 
r: r=100 fjs (*), r=150 /is (*), and r= 200 fxs (t). The 
colors are: zeroth order correlations(black), second order cor- 
relations (red), fourth order correlations (blue), sixth order 
correlations (green), eighth order correlations (orange), tenth 
order correlations (purple), and are twelfth order correlations 
(gray). The signal intensities have been renormalized to em- 
phasize the common long-time decays after ~ 200 /is. The 
lines shown are intended to be a guide and do not represent 
a fit to the experimental data. 



cay at the same rates after ~ 200 fj,s for adamantane and 
~ 100 /xs for calcium fluoride, despite variations in the 
initial density matrix. It has been propose d 25 ' 27 that in 
the asymptotic long-time regime where the FIDs decay 
the same (see Figs. I3I4[) . all elements of the density ma- 
trix share the same eigenvalues and eigcnmodes of decay. 
Our results clearly show that the decay rates appear to be 
the same and, equivalently, that the relative amplitudes 
remain constant. We show the results of the fits of the 
long-time portions of coherence order n to an exponential 
decay 

S n (t)=A n exp(-T n t) (22) 

in Table IIIII for adamantane and calcium fluoride. In 
what follows we give insight into this long time behavior. 

C. Distribution of multispin correlations and its 
dynamics 

The coherence order of the time evolved density matrix 
along the z-axis is conserved - the density matrix only 
connects states whose total moment's projection onto the 
z-axis differs by h (this is a property of high field, secular 
dipolar evolution which conserves the total moment along 
the applied field). The experimentally obtained data re- 
solves build-up and decay of coherence orders transverse 



to the z-axis, which reflect multi-spin correlations in the 
density matrix. With m and n denoting basis states with 
different amounts of transverse magnetization, the quan- 
tity of interest is 

S = Y, ei(n ~ mH \Pmn(t)\ 2 = £ e in *S n (t). (23) 

m,n n=— oo 

Since unitary time evolution preserves the spectrum of 
the density matrix generally and its purity, Tr[p 2 (£)], in 
particular, there is a convenient sum rule Y, n S n (t) = M 2 , 
where M is the initial magnetization of the sample.— 
Also, quite generally, we expect S n (t) to be peaked near 
n = with its width expanding in time (phcnomcnolog- 
ically, it may be approximated as a discrete Gaussian 
profile, e.g. see Fig. [5]). Quite trivially, then, the long- 
time temporal decay of individual S n (t) is dictated by 
the growth of the width of the entire profile and the sum 
rule. Slightly more formally, we can use the second mo- 
ment 

JV 2 (t) = £n 2 Sn(f) (24) 

n 

as a proxy for the distribution's width, and the center of 
the distribution is expected to decay as ~ l/^/N 2 (t). 

The basic qualitative aspect of the dynamics is the 
redistribution of correlations (in the z-basis) from few 
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TABLE III. The fit parameters To and r 2 for the fit of the zeroth and second order quantum coherence decays with time in 
the long-time regime in adamantane and calcium fluoride. The error shown in the table for each value is determined from the 
standard error of the fit. 



Adamantane 


To (ms- 1 ) 


r 2 (ms- 1 ) 


Calcium Fluoride 


To (ms- 1 ) 


T 2 (msf 1 ) 


FID 

100 us magic echo 
150 /is magic echo 
200 /is magic echo 


13.2 ±0.3 
13.2 ±0.3 
13.0 ±0.4 
13.6 ±0.4 


13.6 ±0.3 
13.2 ±0.5 
13.0 ±0.4 
14.2 ±0.3 


FID 

50 fis magic echo 
100 /xs magic echo 
130 ps magic echo 


35.0 ± 1.3 

36.0 ± 1.7 

34.1 ± 1.7 

34.2 ±2.1 


33.7 ±1.2 
35.1 ± 1.0 
34.4 ± 1.7 
32.9 ± 1.2 




T(us) 

FIG. 9. (Color online) Measured even-order multiple spin cor- 
relations in calcium fluoride as a function of time T in Fig. 
[2] The shapes represent the FID (■) and echoes at various 
r: t= 50 fxa (*), r=100 /is (a), and r=130 fia (*). The 
colors are: zeroth order corrclations(black), second order cor- 
relations (red), fourth order correlations (blue), sixth order 
correlations (green) , eighth order correlations (orange) , tenth 
order correlations (purple), and are twelfth order correlations 
(gray). The signal intensities have been renormalized to em- 
phasize the common long-time decays after ~ 100 /is. The 
lines shown are intended to be a guide and do not represent 
a fit to the experimental data. 



spins to the entire crystal. This can be gleaned from a 
set of coupled linear differential equations governing (Li- 
ouvillian) evolution of multi-spin correlation functions, 
known for quite some time^. Crucially, these equations 
are completely specified by the moments of the NMR 
lineshape^. Unfortunately, these moments are not par- 
ticularly well known beyond a few low orders^, so there 
is no closed form solution of this problem in the dipolar 
case of interest here. 

However, some progress has been made recently, 
by Zobov and Lundi n 20 ! 35 , which we now summarize. 
Firslj^, for an infinite range dipolar kernel a wealth of 
closed form solutions was obtained. For instance, 

S n (t) = In(t 2 )(l + M 2 n 2 /t 2 )e- Aht \ (25) 
where M 2 is the second spectral moment of the NMR 



lineshape. The long-time behavior of this model shows 
convergence, albeit non-exponential, of different coher- 
ence orders, S n ~ 1/t for all n as t -> oo (with N 2 ~ t 2 ). 
While this particular feature is reminiscent of our find- 
ings, the infinite range model and many of its proper- 
ties are not physical although with some modifications 
it may be used, perhaps, to describe intermediate time 
behavior.— 

More recently (and realistically), the same group stud- 
ied the dipolar case^S and established simple exponential 
growth of N 2 (t) ~ exp(2Ci). Our data for adamantane 
and calcium fluoride supports this prediction qualita- 
tively, though this measurement is difficult, at least par- 
tially because of the encoding of finite numbers of spins 
in our experiments. However, our fits for the long-time 
portions for calcium fluoride gives a rate constant C on 
the order of 0.02 ^s _1 in calcium fluoride and 0.008 fis^ 1 
in adamantane. Combining this result with an assump- 
tion of a simple Gaussian profile (empirically justified for 
our data but likely more complicated in reality), we may 
expect a simple exponential decay of S n (t) ~ exp(-r„t), 
with approximately constant T n . which is also consistent 
with the data in Figs. El Hand in Tabic Hill 

In connection with predictions based on microscopic 
chaos^I, we note that the decay rates r and T 2 of the 
zeroth and second order correlations matches the FID 
decay rate g in CaFj^. The same is not true in adaman- 
tane, which could be a result of residual motional aver- 
aging in adamantane that makes it a less than ideal test 
case for rigid lattice NMR. In particular, the agreement 
of the decay rates of various coherence orders is consistent 
with the prediction of Ref. H3, which predicts that the 
entire density matrix shares common decay properties in 
the long-time limit. 



V. SUMMARY AND OUTLOOK 

This work aimed to apply novel multiparticlc tech- 
niques for measuring the quality of time reversal in a 
magic echo. In the process, we have provided a new 
quantitative understanding of realistic magic echoes, but 
also caught a glimpse of the fundamental and universal 
process by which closed many-body systems reach ap- 
parent thcrmalization or pseudo-thermalization. Simple 
athermal correlations (easy to produce and to measure) 
evolve unitarily into ever more delicate and difficult to 
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observe multiparticle objects. While the decay of the for- 
mer is commonly observed and equated with thermaliza- 
tion, one of the accomplished goals of the present study 
was to document the latter process. 

It would be of great interest to compare and contrast 
this physics under different conditions and, especially, 
with increased signal-to-noise, e.g. using cryprobes^. Of 
particular fundamental interest is the influence of the en- 
vironmental degrees of freedom - the general expectation 
is that external sources of decoherence in the spin dynam- 
ics will cut off the growth process. In this case the decay 
of low order correlations will be due to extrinsic noise. 
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